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ABS iii Ga: 


The GFCS MK86 solves the ballistics seule for the 577/54 
gun iteratively utilizing table look-ups. Each ballistics parameter 
table is sub-divided two-dimensionally, corresponding to selected 
elevation and slant range intervals. The final value for each para- 
mreter is found by interpolating between parameter values evaluated 
at elevation angles which bracket the actual elevation angle. 

This thesis déseripes the method of determining the ballistics 
solution for AA cartesian ballistics. A description of the ballistics 
geometry and of the model utilized is presented, together witha 
FORTRAN simulation of the MK86 computer program and a listing 


of the ballistics tables used in this system. 
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SARE OMB NCRE Aa 


A. GENERAL DES CR rey 

The GFCS MK86 is the first digital gun-fire control system 
mene considered by the U. 5. Navy. It is designed te conirol the 
5''/54 gun-mount and to solve the fire control problem for the air, 
mummiace, shore-direct or shore-indirect modes of operation, Imihis 


Mayesis Only the AA cartesian ballistics mode is considered. 


B. PROB SEVEDESCRIPEION 

PecOmielter is iimuiattOn proplrain was fell 10 be a prerequisite 
to a complete study and analysis of the GFCS MK86 AA model This 
in turn necessitated a careful study of the existing documentation 
and machine-language software of the MK86 system. 

The FORTRAN simulation of the MK 86 AA mode was divided 
tase three sections: The pre-filter processing section, programmed 
Beer J. D. GORMAN: The three-dimensional numerical filter, by 
LT T. J. WILL; and the ballistics section contained in this thesis. 
These three programs can process radar target information from 
target detection to ballistics output prior to the generation of gun 
orders. 

Another objective of this project was the generation of adequate 
documentation of the systern. This phase was also divided into the 


three sections described in the previeus pagac.co 





er SUMMARY OF INTRODUCTION 

This thesis contains a description of the MK86 ballistics program 
ancluding the ballistics intercept geometry, ballistics medeliceq ang 
sequence of ballistics computations. A description of each subprogram 
indicating its function is also presented. The FORTRAN simulation 
program including the ballistics tables is listed and differences 
between the machine-language program and the FORTRAN program 


are indicated. 








II, INTRODUCTION TO MkK86 BALLISTICS 


The GFCS MK86 ballistics program accepts stabilized and 
fillered target position and rate data, together with ownship's data, 
and computes the gun ballistics required for the generation of gun 


pointing and fuze orders. 


a. mee LbLISTTeo emt Ue 41-2 INITIALIZATION) 

The ballistucsmparanacternsesuch asvaur temapenaturegwair density, 
wind, and initial velocity are entered at the console and @GmPes ponding 
increments from noninal values are computed. 

Target and ownship's data are used to compute the target's 
present position and rate relative to the gun's position at ballistics 
tome. Thiseposition ofwibe target xs calledwepoint 1 (see Figume 1), 
auc! 1SecOmputed inespherical coordinates. The relative, rates are 
resolved into range and crossrange components. 

A test is made to determine whether the input data is associated 
with anew target. Ifilis, an acquisilion mode is declared and the 
processing sequence calls the acquisition subprogram. In the case 
of an old target, track mode is designated and the acquisition sub- 
mropram is bypassed. This is illustrated in the flow chart of 


Mmgure <¢. 
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B. BALLIS TICS ACQUISITION 

Before a ballistics solution can be developed, a first estimate 
of future target position is required. This is called the advance 
Position, or point 3 (see Figure 1). If the system has been tracking 
the target, then the last ballistics solution, updated, serves as the 
advance position. Witha new target, however, some other means 
fer estimating a reasonably accurate point 3 is necessary. This is 
the function of the ballistics acquisition subprogram (see Figure 2). 

im the acquisition pase, point 3 1s computed as follows, “Ihe 
meojectile time of flight to the present target position (point) is 
computed from the ballistics table in the normal way (see Section III). 
Meine this time of flight and the target's present position and velocity, 
an estimate is made of the target's future position (point 3). A number 
of logical tests are employed to check the validity of the estimate. 
Mme slant range, bearing, and elevation of point 3 are designated 


feo, DOYS, and E3 respectively. 


C. DETERMINATION OF POINT 2 

The ballistics program next us es the advance position of the 
target (point 3) as the aim point and computes a predicted "impact" 
point from the ballistics tables. This is called point 2 and the 
corresponding slant range, bearing, elevation, and time of flight are 


designated R2, BY2Z, E2, and T2, respectively. These values differ 


iba 





from the point 3 quantities due to non-nominal powder and air tem- 


peratures, and because of comolhsjedriit, sands) 7 name seca 


im). CAILGULA TION OFS POIeiSs2 ee a7 c- 

The target position which corresponds to the computed timevor 
maent 2 is calculated and designated as point 21 (see Figure 1). Then 
tae differences in slant range, bearing and elevation from those at 
Meme 2 (the predicted projectile impact point) are computed. These 
maerements are designated IR2l1, IBY2l, and lE¢l, respectively, and 
a corresponding linearized estimate of the time of flight increment 
Pr2 is comptted. 

iameliy the pomt clstargem position is updated* by the time 
M_~erement 112 to determine point 22 (see Figure 1), which is the 
iggest EStimate of the point where the target and projectile will coincide 


fameein position and time. 


iB. Pie lo thes ITERATION (1. E. SECOND PASS} 

After point 22 has been calculated once, as indicated in the 
ee eding Sections, the program uses point 42 {(maediiived) as a new 
(improved) point 3 and goes through the calculation to point 22a 
second time. This is shown in Figure ¢. Finally the increments 
Mee2, IBY22, [E22, and IT22 between point 2 and point 22 are 


eomiputed for use in the ballistics output preeeam: 


lez 





Ee. CLOSURE Tisai. 

A test is made at the end of each of the two passes described 
above to determine if the nen point 22 is a valid (i.e. accurate) 
solution. This is done by comparing point 21 data and point 2 data, 
which is called the closure test. The test after the first pass merely 
sets a flag, which is used to modify the criterion used in the test 
which follows the second pass. If the second-pass solution fails to 


close, then the computation aborts for this time interval. 


OF BAULISIIGaeuTeurT 

If closure has occurred, the gun elevation, bearing and time of 
maent which were calculated with point 2 as the aim point are modified 
by the amount necessary to move the impact point to point 22. This 


data is then transmitted to the gun and fuze order programs. 


i 
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lll. BALLISTICS EPELCIS Mepis. 


The MK86 ballistics modeling program provides the means for 
determining the values and effects of the various ballistic parameters. 
mine approach taken and the structure of the modeling is descrabed 
here for AA cartesian ballistics (see Figure 3). The basic function 
of the model is to provide projectile impact or intercept point parameters 
based on advance target position parameters (that is, to provide point 2 
output data from point 3 input data). The ballistics parameters are 


Pbeomen down into the three categories described below. 


A. Le GOmy ONE PARA TERS 
Ni Gun elevation angle kG 


Ns Time of thr nt, be 


a Range €ifect per unmeemmetal velocity variation, DRDU 

4. Range effect perMWmitvair Censity variation, DiKDo 

i. Range effect per unit air temperature variation, DRDTIF 

6. ime .o: went Cifect per whit cneM Vcleeity variation, 
D'TSDU 

ig Time ol might eitect permit ainsdencitavariation, Wis Dr 

oe Dintetom tient ciiecy per wml air volar otmne Varta t1On, 
IBYAS sR S 


ng 





iach of these parameters is modeled by a cubic polynominal 
function of range of the form: 
F=K0+K1(IR3)+K2(IR3)° +K3(IR3)> 
=KO+TR3(K It IRS is 2 eee) 
mevere | is the parameter to be dctermined andwhe K's are coeiiteients 
obtained from tables which are arranged in elevation "stacks", each 
of which is "partitioned" according to slant range. This is illustrated 
by the table below: 


Tig) tye ah 
Typical Ballistics Table for One 


Parameter (suchas Time of Flight) 


Stack. 1 Stack 2 Stack N 
Sitack Stack Same format 
elevation elevation as other 
angle #1 angle #2 stacks 
List of bast of 
range range 
intervals intervals 
in this in this 
stack stack 
Coefficients Coefficients 
Poretirst for -tirst 
range interval range interval 
Coefficients Cocificients 
for 2nd range for 2nd range 
interval interval 
etc. etc. 


The quantity IR 3 in the above cubic polynominal is the incremental 
slant range (i.e. the slant range of interest minus the lower range 


bound of the range partition in use). 


20 





That is, each category-one parameters is modeled by a cubic 
equation in each range-elevation partition, as shown in Figure 4. It 
Heetoe coefficients of the cubics hon are provided by the ballistic 
tables. Inthe simplest case, two values of F are computed, one 
corresponding to each of the elevation stacks which bracket the actual 


elevation. The final value of F then follows from linear interpolation. 


.. CATEGORY. TWO PARAMETERS 

ns Coriolis effects on range, crossrange and time of flight 

oe Drift emkect 

OG Angle of fall factor and calculation of tangent of the angle 

of fall 

Fach of these parameters is modeled by an analytic expression 
modified by a correction factor which is calculated in essentially the 
Same way as a category-one parameter. In this case the factor K is 


eomputed by a linear incremental slant range relation, K=O els). 


C. CATEGORY-THREE PARAMETERS 


Ih lWamee effect per unit preyectile welcnmvariation 
(ae Time Of flight effect per unit projectile weight variation 
oe Wind and motion effects on range and crossrange 


These parameters are calculated in special subroutines which 


do not require range and/or elevation dependent table lookup. 


“ai 





PeOLN 3 
PARAMETERS 


OWN SHIP BE Eis TiC 
EeRA MIST IERS MODEL POINT 2 PAR ee 


B22, Bee 





NON-NOMINAL 
mL UISTICS 
AND WIND __| 
COMPONENTS 


Paetire 3 
BALLISTIC MODEL 


Pee Mie TER 


PARAMETER 
il | HY 


INCREMENTAL INCREMENTAL 
SLANT RANGE SLANT RANGE 


RANGE PARTITION 1 


PARAMETER ‘ 
fl 


INCREMENTAL INCREMENTAL 
SLANT RANGE SANT RANGE 
RANGE PAR TEPOl ie 2 
Bea A Tien BDEVATICN 
STACK #1 SE vCh a2 
Figure 4 


PARAMETER AS A FUNCTION OF ELEVATION AND RANGE 


Ze 





lV QR eee IS TICS en. 


The ballistics tables were derived by modeling tange fae 
data provided LEC by NWL Dahigren [Ref. 1]. There are six sets of 
tables. Each table corresponds to one or more ballistics parameters. 
The AA tables are of the following form: 
i Fach parameter's section of the table is divided into 
elevation stacks. 
rie EFacheelevataenu stack is divided into rangespantitiens 
or intervals. 
S Por cach vampe panmition there are a siet wi coefficients 


which are used to evaluate the chosen parameter. 


Ie MK86 MACHINE LANGUAGE TABLES 

iMticetables are located in the computer im sequential order. | Pie 
address of the first entry for each parameter is stored inthe ETABAD 
matrix. This address is also stored in ISRCH(1). To find the address 
@ietme first elevation stack to be used, the address stored in JISRCH(1) 
fminoditied by the contents of the PUSHTB matrix. The first 3 bits 
@ieecach PUSHTSB entry contains the number ofsentries to be bypassed 
in order to locate the address of the initial elevation stack for the 
low elevation AA mode (advance elevation less than 25). The last 
8 bits indicate the number of entries to be bypassed for the high ele- 


vation AA mode. —The modified address is tmenm stored 11 lon Cue. 


Ce 





The first entry of each clevation stack contains: 


hs Stack elevation angle 
ee, Number or mange partitions, minus one 
om Internal shift count 


The stack elevation angle is stored in STACK(3), the number of range 
partitions is stored inthe NP matrix, and the internal shift count 
Necessary to position the binary point in the KSH matrix. 


A complete listing of the tables may be found in reference 2. 


oe BALLISTICS TABLES IN THE FORTRAN SIMULA TIGN 

The entry into the the organization of the tables is different in 
the FORTRAN simulation program from that in the machine language 
program. 

The ballistics tables are stored in a three-dimensional array. 
The coordinates of one matrix are parameter, clevation stack angle, 
and the location of a particular item in the stack. The first entry in 
a particular elevation stack is the elevation stack angle, and it is 
followed by the number ey range paeriions, Lente limits of the 
range partitions are listed, followed by the coefficients for each 
heamee partition. The location in a stack of the first cocificient 
mequired is calculated by multiplying the number of coefficients times 
the number of the range partition of interest and adding the location 


of the first range partition plus 3. Then the location is stored in T17. 


A listing of the ballistics tables starts on page 90. In this listing, 


24 





each table corresponds to a ballistic parameter (see Table 2), and 


each column corresponds to an elevation stack. 


Table 


10 


lll 


TAREE 2 


Ballistic Tables (FORTRAN) 


Pave tere. 

Timmerotatien, (13) 
ouper-elevation angle (EGB) 
Range effect due to change in 
phitial velocities U ) 

Time of flaght effects dime*to 
Cmanpie i) initial velocity (mo howe 
Range effect due to change in air 
Gensity (DR UH) 

Time of flight effects due to 
change in air density (DT3DH) 
Phomesewcitoet are to change 1n 
Alrabenperaliures( Wil Ul b) 

Time of flight effects due to change 
inal r"temperatTre (Die TE) 
eangemt ofin cle wr fall (TANW ) 
Coriolis effects 


Drift effects 





V. BALLISTICS SUBPROGEA i> 


Since this project was concerned only with the AA aspect of the 
GFCS MK86 system, all references to surface ballistics, grid mode, 
ROS mode, TDT mode, and spherical ballistics have been omitted. 
imereine flow chart of the subprograms see Figure 5. A complete 


explanation and flow chart may be found in References 3 and 4. 


a, die Gl 
imis subpregcram ireezes or holds the console inputs “ie the 


ie 


duration of the % second ballistics computation cycle, and calculates 


and sets the following indicators: 


iL: AA ballistics 

ay Cartesian ballistics 

on Aimimno type, charge 

4. Gun and gun control console assignments 
5. Target assignment 

op Pieqmisition or tracking ballistics 


in addition, LAIG] coniputes mon neminaisinereimments thatginay 
exist in various ballistics parameters, performs a ship's course 
update, and determines gun parallax corrections. Calculations are 
also performed to update the target data to “ballistic time” (1. €., the 


next "integer''Z second). 
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.. eee 

This subprogram performs the acquisition ballistics computations 
for AA cartesiar ballistics. Its function 1s to calculate point > (the 
theoretical aiming point) from point | (the target position at ballistics 
time). This subprogram determines the first advance position (aim 
point) when no valid prior advance position exists. This situation 
fould occur whenever: 

i A new AA target has been assigned 

Gi There was no ballistics solution during the previous 


4-second ballistics computation cycle 


Bi There has been a change in the mode assignment for 
the son 
4. There has been a change in the ammo (ype cigar c 


The countdown flag (KGOOD) is set to one, indicating that the 
meamisilion subprogram must be completed more than once before the 
eoimers released from any previous state it may have beenin. The 
meaquire flag is set is one to indicate {he acquisition ballistics miode, 

The advanced elevation (E3) 1s set equal to the elevation angle 
at point 1 (E11). TMONIT, which controls the ballistic table computa- 
tions, is then entered to find the nominal time of flight (T!) from RC, 


ime ranpe to point f. 


ca 





C. BCOMPA 

This subprogram computes the time of flight to pois (73) 
and its derivatives with respect to range and elevation. It also 
Semiputes the mixed second partial derivative with respect tovramgec 
and elevation. Using E3 and R3 as inputs, TMONIT is called and 
the time of flight is returned in TEMP2, the range derivative (DT3DR) 
in TMP3, the elevation es (DESDE) in DPDEP, and the miccaa 


partial (DT3RDE) in DPDED. 


iD BCOMPB 

This subprogram computes the super elevation angle (EGB) and 
its derivatives with respect to range and elevation, as well as the 
mixed second partial with respect to range and elevation. TMONIT 
is called and the super elevation angle is returned in TEMPZ, the 
memee derivative in TMP3, the elevation derivative (DEGDE) in 
ieee, and the mixed partial (DEGRDE) in DPDED. Then EGB and 
E3 are summed to form the total elevation angle E4=E3+EGB, which 


feemmuted to a maximum of 85 degrees. 


dig BCOMPC 

This subprogram computes the derivative of range with respect 
fora nOn-nominal increment of projectile velocity (DRDU) and the 
derivative of DRDU with respect to elevation (DRDUDE). Then the 
metal tame Glace, (URI die te nora oneal projectile velocifyeis 


COMaputed. 
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i. iG@ Mig 

This subprogram computes the derivative of {he fimic ot aan 
math respect to the increment in projectile velocity (DR3w wae ne 
meese of flight effect (UT3) due to nen-nominal projcciile 1a - 
The mixed second partials of time of flight with respect to increment 
mmerojectile velocity and elevation (DT3UDE) and range (DT3U DR) 
are also computed. 


~~ 


CP BCOMPE 

This subprogram computes the derivative of range with respect 
to non-nominal increments in air density (DRDH) and the derivative 
oe ORMDriewith respect to elevation (DRDHDE). Then the actual range 


ermeet due to non-nominal ballistic air density (HRH) is computed. 


ie BeOMPr 

This subprogram computes the derivative of time of flight with 
metepect to a non-nominal increment in air density (DT3DH) and the 
actual time of flight effect (HT3) due to DT3DH. The mixed second 
pee iials of fume of flight with respect to air density and range and 


efevanon (DT3HIDR, DI3HDE) are also computed. 


1b BCOMEG 
This subprogram computes the derivative of the range with 


Pespect to a non-nominal increment in surface air temperature (Ol), 
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the actual range effect (TFR3), and the mixed partial of range i 


respect to temperature and elevation (DRTFD£). 


ie BCOMPH 

This subprogram computes the derivative of time of flight with 
feeeeect tO a non-nominal increment in surface temperature (D?P3DI 8); 
Semeputes the actual time of flight effect (TFT3) due to air temperature; 
and calls the weight effect subprogram (RTWT). 

The subprogram then performs a partial sum of all the range 
effects and time of flight effects calculated so far and calculates the 
derivative of the time of flight to point 2 (the predicted impact point) 


with respect to range (DTZDR). 


K. me COV | 

This subprogram computes the downrange and crossrange 
components of gun motion (DMR30G DMGOG) before calling GUNMOT, 
freweun motion subprogram. It then computes the tangent of the angle 


of fall (TAN W) and calls the wind effects subprogram, WIND. 


i. BCOMPI 
This subprogram obtains a first partial sum of range effects, a 
second partial sum of time of flight effects, and calls CORLISS to 


obtain coriolis effects. 
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M. nCOMIPK 

This subprogram computes the intercept geometry ty Gereriii ie 
the horizontal range, vertical range, slant range, and bearing to the 
predicted projectile impact point (point 2). It also computes the 


elevation angles for the gun. 


N. TARG3 

This subprogram performs the iteration and closure computations. 
Its function is to update through two iterations the previous advanced 
aim point (point 3) such that satisfactory coincidence is obtained 
between the predicted impact point (point 2) and the predicted target 
position at impact time (point 21). The closure equations determine 
meren there is sufficient closeness between the two points for an 
acceptable solution. 

The input ice to [AR G3 come from two sources, emrnic 
tracking mode, the data from the previous pass through TARG3 to 
mBCOMIPK is used to update current values. In the acquisition mode, 
the values calculated in a eee 22a throu heb GO (ame 


used. 


oF TARG4 
This subprogram computes the ballistics outputs. It computes 
the final inertial gun orders, gun order rates, and fuse time rate as 


well as the inertial impact coordinates and rates for eventual display 


oy 





on the B-scan. These computations are oriented to the start of the 


next 1/16 second time-slice. 


le TMONIT 

For a given parameter, this subprogram does the initial setup 
and locates the proper table and elevation stack corresponding to the 
target's elevation. Then it calls appropriate subprograms to cal- 
eilate the value of the parameter from the chosen table. 

In more detail, either two or three stacks are selected so that 
Beir elévation values bracket the value of E3. Two are selected for 
parameters which are linearly interpolated with respect to elevation 
and three for those which are parabolicaly interpolated. The partition 
selection, coefficient look-up, and polynominal evaluation with respect 
pomtange are performed on each of the selected sets of coefficients 
yielding a value of the parameter and its derivatives corresponding 
to each of the bracketing elevations. The final value of the para- 
meter and its derivatives are then obtained by interpolating with 


respect to elevation. 


©. AASRCH 

this subproeraim controis’ the selection of coectiicicnrsmion 
bracketing elevation stacks, the cubic polynominal evaluation, and 
the interpolation between elevation stacks. It also controls the 


calculation of the appropriate derivatives. For standard tables (four 


Be 





coefficients per partition) the address vf ithe first coefficient is found 

by calling the SEARCH subprogram. Then the polynomial subprogram 
temcalled to evaluate the parameter, which is stored, with its derivative, 
in the STACK matrix. When all parameters and derivatives correspond- 
ing to the selected elevation stacks have been evaluated, then the inter-~- 
polation subprogram is called to get the final values. For non-standard 
tables, only the address of the first coefficient is found for each selected 


Pievation stack. IJhese values are stored in the HOLDB3 matrix. 


rs oRARCH 
Minis subprocram, which is called by AASR CH finds ihe correc: 
range partition and the location of the first coefficient for that range 


partition. 


S. POLYN 
This subprogram, which is called by AASRCH, evaluates the 


polynomial from the coefficients and finds the derivatives with 


mespect to range for each elevation stack chosen. 


A. eR 

This subprogram performs linear or parabolic interpolation on 
a set of two or three parameter values with respect to elevation. 
the interpolates value is the value of the parameter at the actual 


advance clevation. 





i... DECSti 
This subprogram decodes the elevation information from each 


mack and stores il ino l2CMia). 


i. SLIDE 
Thisssubprogranmmmoves the paxamctersystored ingihess 0.0m 


Beatrix up one location. 


W. FINDAL 
This subprogram finds the location of the next available elevation 


stack. 


De CORLIS 

This subprogram computes ballistic Coriolis effects for slant 
range, allitude, time of flight, crossrange and downrange effects. 
The Coriolis table is a non-standard table with six coefficients per 
femese partition. Three linea coefficients are computed and used to 


Sreterimine the effects of Coriolis. 


. DRIFT 

This subprogram computes the crossrange effect of projectile 
mm. Lhe drift effect table is a non-standard table with two co- 
efficients per range partition. Two polynomials are evaluated for 
AA ballistics, one for each elevation stack bracketing the actual 


advance elevation. 
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Zi GUNMOT 

This subprogram calculates the ballistic effect due to gun 
motion (i.e. the advance range, time of flight, deflection and cross-— 
range effects). Under conditions of no gun motion, the projectile 
leaves the muzzle with an initial velocity vector (U) at inertial true 
Bearing BCY, and inertial elevation, £G. With gun motion, the 


Somaponent along the U vector will increase the magnitude of U. 


ee RIWT 
This subprogram calculates the (minor) effect ol projectile 

weight on time of flight and range. Due to inertia, an increase of 
projectile weight will produce a decrease of initial velocity. It will 
also produce a decrease in drag acceleration, which is equivalent to 
a decrease in air density. Thus the range and time of flight effects 
due to a projectile weight variation can be computed as the sum of 
the equivalent effects of increments in initial velocity and ballistics 


paliaensity. 


AB. WIND 
This subprogram calculates the total effects of wind on advance 


range, time of flight, and crossrange. 


fee. ALIMIT 
This subprogram limits the Z coordinate ofa target toa 
minimum value, based upon target horizontal range, earth's curvature, 


and height of the gun above the sea surface. 
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VI. FORTRAN SIMULATION 


Since only the AA mode was to be considered, all ilags and 
variables necessary to indicate the AA mode were pre-set to insure 


that only the AA cartesian ballistics computations were performed. 


ae ASSUMPTIONS MADE 
The FORTRAN program has no real time clock, therefore 
certain timing assumptions had to be made. These assumptions are: 
a JOELCK (target valid data-time for ballistics) 


equals -1/32 second 


Cc. TBR (saved value of ballistic time) was initialized 
to 0 
3. TB (ballistics time to next second) initialized to 


% second 
4. TBFIN (final ballistics time for next 1/16 second 


Vaimadlimec te 15/82 eeecme 


loa FORTRAN PROGRAM 
The program, which is listed starting on page closely follows 
the arrangement of the MK86 machine language program. It has 
been suitably annotated with comment cards, and otherwise should 
be self-explanatory. It has been written so that all ranges are in 


yards, and all angles are in degrees. 


4] 





Three internal subprograms (DI¢CSI, SLIDE, AND FINDAL) 
of ITMONIT in the machine language program have been treated as 
separate subprograms inthe FORTRAN program. Also, in BCOMPB 
the super-elevation angle is computed in minutes and then converted 


to degrees. 


cS. VERIFICATION 

Complete verification of the FORTRAN program was not done 
due to lack of sufficient data. The separate subprograms were 
verified by calculations done by hand. Using data furnished by LEC, 
pepertial verification was attempted. However, effects such as wind, 
coriolis, non-nominal parameters and ship's motion had been set to 


zero in the data provided [Ref. 5]. 
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